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In this note we consider three issues related to the unitary Fermi gas in a harmonic 
trap. We present a short proof of a virial theorem, which states that the average 
energy of a particle system at unitarity in a harmonic trap is twice larger than the 
average potential energy. The theorem is valid for all systems with no intrinsic scale, 
at zero or finite temperature. We discuss the odd-even splitting in a unitarity Fermi 
. gas in a harmonic trap. We show that at large number of particles N the odd-even 

splitting is proportional to N l < 9 huj, with an undetermined numerical constant. We 
also show that for odd N the lowest excitation energies are of order N^^ftw. 

Recently Fermi gas at unitarity has been realized and studied experimentally. In most 
experimental setups the system is confined in a harmonic trap, hence of great interest are 
the properties of such a system in such traps. In this note we explore several aspects of this 
system: the virial theorem, the odd-even splitting, and the excited energy levels of a system 

Ch ' with an odd number of particles. 

O ' 

1. Reference [1] contains a proof of a virial theorem for a unitary Fermi gas in a harmomic 
trap. This theorem states that the total energy is exactly twice the average potential energy 
due to the trap, (H) = 2(V). The proof is based on the force balance in a harmonic trap 
l/"") ! within the local density approximation. There is, however, a simpler proof based on the 
Hellmann-Feynman theorem. It is similar to, but arguably simpler, than the a proof due to 



^ . Chevy quoted in Ref. [2]. 
O 

O 



Consider the spherical trap first. The Hamiltonian of the system can be written as 



N 



a=l 



2m 



X. 

where i/2-body includes all two-body interactions. We have highlighted the fact that the 
Hamiltonian depends on the trap frequency uj. 

Since in the unitarity limit the system lacks an intrinsic scale, the ground state energy 
has to be proportional to fko, which is the only energy scale available. Denote the ground 
state of N particles as \u, N), we have 

(u,N\H(u)\u,N) = c N fitJ, (2) 

where cn is a constant dependent on N. 

The Hellmann-Feynman theorem states that 

u-^(u,N\H{u)\l;,N) = (u,N\u^^\uj,N). (3) 
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We see from Eq. (2) that the left-hand side is equal to the ground state energy, and from 
Eq. (1) that the right-hand size is equal to twice the average potential energy. One obtains 
the virial theorem: (H) = 2{V). 

It is easy to extend the theorem to the case of an anisotropic trap. In this case the 
Hamiltonian is 



N 3 
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" = ££(S + ^)+'w. <*> 

a=l i=l 

and instead of Eq. (2) we have 



{u i ,N\H(u i )\ui,N) = E N (u i ). (5) 
Consider a generalization of Eq. (3): 

Y,^(^,N\H(u t )\u t ,N) = ( Ui ,N\^uJ¥^\ui,N). (6) 

i i 

Basic dimensionality analysis tells us that E^iui) is equal to a common frequency, e.g., the 
geometric mean (cji^^) 1 / 3 , times a dimensionless function of the ratios between frequen- 
cies. Such a function is homogeneous function of first order of its arguments, i.e., 

3 d 

^2u)i—E N (u)i) = E N (ui). (7) 
i=i 1 

while the Hamiltonian itself is a homogeneous function of second order with respect to lo^. 
From Eq. (6) we obtain the virial theorem in an anharmonic trap. 

The proof can be extended to finite temperature T as well. For simplicity consider the 
harmonic trap. The free energy has to scale as 

F N (u,T) = c N (^Jhcu. (8) 

Now can depend on the dimensionless ratio T/uj. We now apply the finite-temperature 
version of the Hellmann-Feynman formula, 

f) r)M 

w _^,r, = <^>. (9) 

The right hand size is again 2{V). The left hand side, from Eq. (8), can be transformed to 

which, according to a thermodynamic relation, the average energy (if) in the canonical 
ensemble. The extension to anharmonic traps and grand canonical ensemble is straightfor- 
ward. 

We note that the proof does not rely on the local density approximation. Compared to 
Chevy's proof (quoted in Ref. [2]), here we rescale the potential while in Chevy's proof the 
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wavefunction is rescaled. The theorem holds for systems with spin imbalanced, independent 
of their phase structure. 

Finally we note that the virial theorem should be valid for other systems with no intrinsic 
scales, for example for the system of anyons in two-dimensional harmonic traps. 

2. The quantum Monte-Carlo calculation of Ref. [3] shows a clear even-odd effect, rem- 
iniscent of the behavior of the nuclear binding energy. If one denotes by En the ground 
energy of N particles in an isotropic harmonic trap (E N = c^hw) then the quantity 

Aat = E N — -(E N ~i + E N+ i), (11) 

for odd N, is positive and approximately 0.6 — 0.7ftw for the range of particle number N 
studied (3 to 21). Recall that the ground state for even N has equal numbers of spin-up 
and spin-down fermions, while for odd N the numbers differ by one. 

The effect is clearly related to pairing between particles, but a question arises: will the 
odd-even effect remains constant when one goes to large iV? Here we argue that in the limit 
of large N, the even-odd splitting has to grow as a small power of N: 

A N = E N - ^(E N -i + E N+l ) ~ N 1/9 u, N large and odd (12) 

To see how does this dependence arise, first recall that in the local density approximation 
(LDA) the system, with an even number of particles, can be thought of as a Fermi gas with 
a spatially dependent chemical potential: 

fj.(r) = (jlq - V(r) = (jlq - imwV (13) 

(we consider an isotropic trap). Assume we have a system with even N, and discuss the 
process of introducing an extra particle into it. For N >> 1, the extra particle can be thought 
of as an extra fermionic quasiparticle carrying definite spin. If N >> 1 the cloud is large, and 
the quasiparticle can be approximately localized at a position r. The energy cost of doing 
so is the chemical potential, plus the local energy gap A(r) (due to superfluidity), which is 
proportional to jj,: 

A(r) = C/i(r) (14) 

(C ks 1.2 according to a Monte-Carlo simulation [5]). Therefore the extra particle will be 
localized where the gap is smallest, which is near the edge of the cloud, r = R. Unfortunately, 
the LDA breaks down there. 

Since the quasiparticle is localized on the thin shell, the curvature of the shell can be 
completely ignored. As the result, the problem is mapped onto the problem of inserting a 
quasiparticle in a symmetric system at /i — in the linear potential 

V(r)=Sz (15) 

where 



S = muj 2 R. 



(16) 
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FIG. 1: Localization of an extra unpaired quasiparticle on the surface of the cloud. The solid line 
is the effective chemical potential, and the dashed line is the effective potential experienced by an 
extra unpaired quasiparticle. The picture is not quantitatively accurate near the edge of the cloud. 

The choice of the notation reflects the fact that £ can be interpreted as an electric field. 

In the potential (15) the symmetric system, in the local density approximation, fills the 
z < half of space; the other half z > is empty. Beyond LDA there is a boundary layer 
around z — 0. The density tends to infinity as z —>■ — oo but it is fine for a scale-free system 
as the unitary Fermi gas. Assuming £ > 0, the quasiparticle cannot stray toward negative z, 
because the gap there is large, and it cannot stray toward positive z because of the potential. 
Therefore it has to be localized at the boundary z = (Fig. 1). At unitarity the ground 
state energy of this quasiparticle can be estimated by dimensional analysis by constructing 
the unique combination of h, m, and £ with the dimension of energy, 

£ extra = • (17) 

Here \ is a universal dimensionless constant. 

Recall that the size of the cloud R is related to the particle number N by 



i2 = e 1/4 W— (24iV) 1/6 (18) 
V muj 

where £ is an universal number (denned as the ratio of energy of an unitary Fermi gas and 
the energy of a free Fermi gas with the same density), we find, by combining Eqs. (16), (17), 
and (18), 

A N = X £ 1/6 {24N) 1/9 hu (19) 

which is the behavior advertised in Eq. (12). Thus we conclude that the odd-even splitting 
grows as a small power of N at large N. 

While £ has been evaluated by many methods (which typically give £ ~ 0.4 [4]), we have 
no previous evaluation of x- Ideally, one would like to use quantum Monte-Carlo method to 
find it, but it should be also possible to estimate x by using the e expansion technique [6]. 
One could use the Bogoliubov-de Gennes equation to find a rough estimate of x- This is, 
however, beyond the scope of this note. 
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We can also discuss the odd-even splitting in an anharmonic trap. We take the limit of 
large N, fixed aspect ratios of the trap. The extra particle is still located near the edge of 
the cloud, but instead of spreading out over the whole ellipsoidal edge it will be attracted 
to the areas where the potential gradient (the effective electric field £) is smallest. This 
corresponds the locations farthest from the center. Repeating the calculations we shall find 

A N = ^6(24^)1/9^1/3^/3 (2Q) 

where u> = (c^iCc^s) 1//3 and cj min is the minimal among cui, and lu 3 . 

3. Let us now discuss the excited energy levels of a system in a harmonic trap. For 
a system with even N and equal numbers of spin-up and spin-down fermions, the lowest- 
energy excitation can be thought of as a coherent excitation of the superfluid, and the lowest 
energy excitations have energies of order fkv [7]. 

For an odd system, another type of excitations exist. Recall that we now have a extra 
quasiparticle moving on the surface of the cloud. We assume that this quasiparticle has the 
dispersion relation e = k 2 /2m* at small momentum k along the surface, where m* is an 
effective mass. In the ground state the extra particle has momentum zero. One can give 
this quasiparticle a nonzero orbital momentum, as the result one get a series of excitation 
level whose energy depends on the orbital momentum £ as 

= hH{£ + 1) = £(£ +1) 1 m 
1 2m*R 2 2£V2 (24AT)V3 m * ^ > 

What is interesting is that the splitting between lowest energy levels is not given by the 
scale f\w as in even systems, but by a much smaller scale N^^huo. 

It is also interesting to discuss the excitations of an odd system in an anharmonic trap. 
Again we consider the limit of large iV at fixed trap aspect ratios. If the smallest frequency 
is unique (e.g.., if uj\ < 002 and uji < 00%) then the extra particle is attracted to two opposite 
points on the longest axis of the ellipsoidal cloud. When the cloud is large the particle can 
jump between the two points only by quantum tunneling. Thus for odd number of particles 
N the ground state is almost degenerate with the first excited state, which has an opposite 
parity. This "parity doubling" does not hold when the trap is an oblate spheroid (e.g., if 
uj\ = uj 2 < UJ3), since now the extra particle is concentrated to a circular ring around the 
cloud edge. 

I am indebted to G. Bertsch, A. Bulgac, S. Y. Chang, J. E. Drut, M. Forbes, A. Kryjevski, 
and S. Tan for discussions leading to this note, and to J. Thomas for comments on the 
manuscript. This work is supported, in part, by DOE Grant No. DE-FG02-00ER41132. 
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